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". Abstract 

o , 

We investigate several quantities, denned in the decays of top quark pairs, 
which can be used to explore non-standard Wtb interactions. Two new angular 
asymmetries are introduced in the leptonic decay of top (anti)quarks. Both are 
very sensitive to anomalous Wtb couplings, and their measurement allows for a 
precise determination of the W helicity fractions. We also examine other angular 
and energy asymmetries, the W helicity fractions and their ratios, as well as spin 
correlation asymmetries, analysing their dependence on anomalous Wtb couplings 
and identifing the quantities which are most sensitive to them. It is explicitly 
shown that spin correlation asymmetries are less sensitive to new interactions in 
the decay of the top quark; therefore, when combined with the measurement of 
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other observables, they can be used to determine the tt spin correlation even in 
the presence of anomalous Wtb couplings. We finally discuss some asymmetries 
which can be used to test CP violation in tt production and complex phases in 
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the effective Wtb vertex. 

1 Introduction 

Precision studies have been in the past a powerful tool to explore new physics at scales 
not kinematically accessible. With the operation of the Large Hadron Collider (LHC), 
top physics will enter into the era of precise measurements [1]. Due to its large mass, 
close to the electroweak scale, the top quark is believed to offer a unique window 
to physics beyond the Standard Model (SM). New interactions at higher energies may 
manifest themselves in the form of effective couplings of the SM fermions, especially for 
the top quark, much heavier than the rest. In this work we concentrate ourselves on the 
Wtb vertex. Within the SM this coupling is purely left-handed, and its size is given by 
the Cabibbo-Kobayashi-Maskawa matrix element V t b, which can be measured in single 
top production [2-4]. In new physics models, departures from the SM expectation 
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V t b — 1 are possible [5,6], as well as new radiative contributions to the Wtb vertex [7,8]. 
These corrections can be parameterised with the effective operator formalism. The 
most general Wtb vertex containing terms up to dimension five can be written as 

C = --£=b-f(V L P L + V R P R )tW- 

iif {9lPl + 9rPr) 1 w; + hx - ' (1) 

with q = Pt — Pb (we follow the conventions of Ref. [9] with slight simplifications 
in the notation). If CP is conserved in the decay, the couplings can be taken to 
be realQ Within the SM, Vl = V t b — I and Vr, gi, 9r vanish at the tree level, 
while nonzero values are generated at one loop level [10]. Additional contributions 
to Vr, gL, gR are possible in SM extensions, without spoiling the agreement with 
low-energy measurements. The size of a Vr term is constrained by the measured 
rate of Br (6 -> s 7 ) = (3.3 ± 0.4) x 10" 4 [11]. A right-handed coupling \V R \ > 0.04 
would in principle give a too large contribution to this decay [12] which, however, 
might be (partially) cancelled with other new physics contributions. Hence, the bound 
\Vr\ < 0.04 is model dependent and does not substitute a direct measurement of this 
coupling. Similar arguments applied to the a^ v terms do not set relevant constraints 
on gR, because its contribution is suppressed by the ratio q v /Mw for small q u . 

Top production and decay processes at LHC allow us to probe the Wtb vertex [2,4, 
9,13,14]. Top pair production takes place through QCD interactions without involving 
a Wtb coupling. Additionally, it is likely that the top quark almost exclusively decays 
in the channel t — > W + b. Therefore, its cross section for production and decay gg, qq — * 
ti — > W + bW~b is insensitive to the size and structure of the Wtb vertex. However, the 
angular distributions of (anti)top decay products give information about its structure, 
and can then be used to trace non-standard couplings. Angular distributions relating 
top and antitop decay products probe not only the Wtb interactions but also the spin 
correlations among the two quarks produced, and thus may be influenced by new 
production mechanisms as well. On the other hand, single top production is sensitive 
to both the size and structure of the Wtb vertex, involved in the production and the 
decay of the top quark [4, 13, 14]. 

In this paper we explore the sensitivity of several quantities, like angular and en- 
ergy asymmetries, helicity fractions and ratios, to new non-standard Wtb interactions. 



1 A general Wtb vertex also contains terms proportional to {pt + pb)^, and a^^pt +Pb)v Since 
b quarks are on shell, the W bosons decay to light particles (whose masses can be neglected) and the 
top quarks can be approximately assumed on-shell, these extra operators can be rewritten in terms 
of the ones in Eq. (T]) using Gordon identities. 
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Although these observables are theoretically related, the experimental determination is 
more precise for some of them than for others. In particular, the experimental precision 
is dominated by systematics already for a luminosity of 10 fb _1 , and a good choice of 
observables can improve significantly the limits on anomalous Wtb interactions. Our 
analysis here is kept at a purely theoretical level, identifying the quantities which are 
a priori more sensitive to anomalous couplings, and estimating the precision in their 
experimental measurement from a detailed simulation, which has been presented else- 
where [15, 16]. 



2 W helicity fractions and ratios 

The polarisation of the W bosons emitted in the top decay is sensitive to non-standard 
couplings [17]. The W bosons can be produced with positive (right-handed), negative 
(left-handed) or zero helicity, with corresponding partial widths T R , T L , T , being 
r = T(t — > W + b) = T R + T L + T . The T R component vanishes in the mj = 
limit because the b quarks produced in top decays have left-handed chirality, and for 
vanishing m b the helicity and the chirality states coincide. The three partial widths 
can be calculated for a general Wtb vertex as parameterised in Eq. (DQ), yielding 

2 1 -*| ( 2 

+ [\9l\ 2 + \9r\ 2 ] (l-Xyy + xl) -Ax b Reg L g* R 
-2-^ Re [V L g* R + V R g L ) (l - x 2 w - x 2 b ) 



+2 lt Xb Re [Vl9 * L + VR9 * R] ^ + ^ " ^ } ' 
2 1 -* I 

\l = ^{[\V L \ 2 + \V R \ 2 ](l-x 2 w +4)-4x h RsV L V 1 



m 2 



Ofi'il 2 + \9r\ 2 ] (l-x 2 w - 2x 2 b ~ x 2 w x 2 b +xt) - 4x b Reg L g* R 



M 2 W 

-2^- Re [V L g* R + V R g* L ] (l - x 2 w - x 2 b ) 
+2-^-x b Re [V L gl + V R g* R ] (l + x 2 w - x 2 ) } 

^4 m " lVRl2] + °" L ' 2 " (1 " 

+2x w Re [V L g* R - V R g* L ] + 2x w x b Re [V L g* L - V R g* R }} 

X (l 2X^T 2Xfo + Xyy 2XyyX[ J + x b ) , (2) 
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being X\y = M w jm tl x b = m b /m t and 

I? I = T-(™>t + M w + < ~ 1rn 2 t M^ - 2m 2 t m 2 b - 2M 2 v m 2 b ) 1 ' 2 (3) 

the modulus of the W boson three-momentum in the top quark rest frame. The total 
top width is 

21 -*| 2 

g \q\ m t rn |2 , |21 / 2 r> 2 4 , _2 ji 



32tt 



{ [|Vi| 2 + |Vr| 2 ] (1 + x 2 w - 2x 2 b ~ tew + AvA + 4) 



-12x 2 w x b Re V L V* R +2 [\g L \ 2 + \g R \ 2 ] (l - ^ - 2x 2 b - ± - ^ + x*) 

-\2x 2 w x h Reg L g* R - 6x w Re [F L ^ + V R g* L ] (l - - xf) 

+6x w x b Re [V L ^2 + V R g* R ] (l + x 2 w - x 2 b ) } . (4) 

The different polarisation states of the W boson are reflected in the angular distribution 
of its decay products. Let us denote by 6} the angle between the charged lepton three- 
momentum in the W rest frame and the W momentum in the t rest frame. The 
normalised differential decay rate for unpolarised top quarks can be written as 

' <!l l{l + cos6}) 2 F R + l{l-cos6*z) 2 F L + ~sm 2 e* t F , (5) 



T d cos 9* t 



with Fi = Ti/T the helicity fractions. The three terms correspond to the three helicity 
states, and the interference terms vanish [18]. At the tree level, F = 0.703, Fl = 
0.297, F R = 3.6 x 10~ 4 , for m t = 175 GeV, M w = 80.39 GeV, m b = 4.8 GeV. The 
resulting distribution is shown in Fig. [U calculated from the analytical expressions in 
Eqs. ((21)-© and also with a Monte Carlo simulation. The latter is performed using 
our own ti generator, which uses the full resonant matrix element for gg, qq — > ti —>■ 
W + bW~b — > fifibfif^b, and hence takes into account the top and W widths, as well as 
their polarisations. Anomalous couplings in the decay may also included in the event 
generation. We observe that finite width corrections have a negligible influence in the 
distribution, and hence Eqs. ((2J) _ © can be used to make precise predictions for the 
distributions. 

The good agreement between the analytical calculation (with the top quark and 
W boson on their mass shell) and the numerical one can be explained substituting 
m t -> m t (l + (tT-t/mt), M w -> M w (l + £ W T W /M W ), with £ t , £ w of order unity, in 
the expressions of the helicity fractions (here we introduce subscripts to distinguish the 
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Figure 1: Differential distribution in Eq. ((5|) within the SM, calculated analytically and 
with a Monte Carlo simulation. 



top quark and W boson widths). We obtain 
F = 0.703 



I (,.097 ( 6— ) - 0.595 Uw^f- J - °' 545 ( & — 
m t J V M w 1 V m t 



+0.055 + 0.487 U£L)(, w ^ l+ ... 



Fn = 3.6 x 10" 4 



I i. is ( &I±) + 2.48 f^^l + 13-22 fe*— 
mtJ V AfW V m * 



+ , 77U „^)- 12 , 8 (^)(^, + ... 



F L = 0.297 



I 1.107(6— )+ 1-405 f^^ 1 
m t J V m h/ 



1.27! ; 6 — 



-0,34,^) _ 1 ,37(^)(^, + 



(6) 



Linear terms have no effect when integrated with symmetric Breit-Wigner distributions, 
and the quadratic terms are very small. 

In the presence of anomalous couplings, the helicity fractions Fj are modified with 
respect to their SM values quoted above. Their variation is plotted in Fig. [2], considering 
that only one coupling is different from zero at a time and restricting ourselves to the 
CP-conserving case of real Vr, qr and g^. We observe that F L and F are much more 
sensitive to than to g^ and Vr. This is due to the interference term V^g^ , which 
is not suppressed by the bottom quark mass as for the gi and Vr couplings. This 
linear term dominates over the quadratic one and makes the helicity fractions (and 
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Figure 2: Dependence of the helicity fractions Fi = Ti/T on the anomalous couplings 
in Eq. (jTJ), in the CP-conserving case. 



related quantities) very sensitive to qr. We also remark that the phases of anomalous 
couplings influence the helicity fractions through the interference terms which depend 
on the real part of Vr, gi and gR (we have taken Vl real, and normalised to unity). 
Thus, the effect of complex phases is specially relevant for where the interference 
term dominates. In any case, the maximum and minimum deviations on the helicity 
fractions are found for real, positive and negative (not necessarily in this order) values 
of Vr, gR and g^. The possibility of complex couplings is examined with more detail 
in section El 

The helicity fractions can be experimentally extracted from a fit to the cos Q\ distri- 
bution using Eq. j5]). In order to estimate the limits on anomalous couplings that can be 
set from their measurement, we assume that the central values obtained correspond to 
the SM prediction, and take their errors from Refs. [15,16,19], giving F ~ 0.703±0.016, 
F R ~ 3.6 x 10" 4 ± 0.0045, F L ~ 0.297 ± 0.016. For these values, it is found that F 
and F L have a similar sensitivity to gR, while the dependence of F R on this coupling 
is smaller. On the other hand, the measurement of Fr sets the strongest constraint on 
Vr and gL- The resulting bounds are summarised in the first column of Table [U These 
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Figure 3: Dependence of the helicity ratios pr ; l = ^r,l/^q on the anomalous couplings 
in Eq. (DQ), in the CP-conserving case. 

and the rest of limits throughout this paper have been obtained with a Monte Carlo 
method, as described in appendix [El 

F\ Pi 

V R [-0.062,0.13] [-0.029,0.099] 
g L [-0.060,0.028] [-0.046,0.013] 
g R [-0.023,0.021] [-0.025,0.026] 

Table 1: la bounds of anomalous couplings obtained from the measurement of helicity 
fractions Fj and ratios pj. 

The sensitivity achieved for non-standard couplings may be greater if we consider 
instead the helicity ratios pr } l = ^r,l/^o = Fr^/Fq, shown in Fig. [3]§| These ratios 
can be directly measured with a fit to the cos 6^ distribution as well. From the expected 
precision in their determination in Ref. [15, 16], and assuming that the central values 
correspond to the SM prediction, we have p R ~ 0.0005 ± 0.0026, p L ~ 0.423 ± 0.036. 
From these values, the limits given in the second column of Table CD can be obtained, 
with an important improvement for Vr and g^. As it has been remarked in the in- 
troduction, the reason for the improvement is that systematic errors, which dominate 
the precision of the measurements (see Refs. [15,16] for details), are much smaller for 
helicity ratios than for helicity fractions. 

To conclude this section, we would like to stress the importance of keeping the 
bottom quark mass in the calculations. Within the SM the nib correction to the helicity 

2 We note that, for a better comparison among them and with other observables, the scale of the 
y axis in each plot is chosen so that the range approximately corresponds to two standard deviations 
(with the expected LHC precision) around the theoretical SM value. 
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fractions is small, of order x\ = 7.5 x 10~ 4 , as it can be seen in Eqs. ([2]). However, as it 
can also be observed, the interference terms involving gi or Vr couplings with Vl are 
proportional to Xb = 0.027, and are of similar magnitude as the quadratic terms. The 
effect of including nib in the computations is illustrated with more detail in appendix 
lAl Nevertheless, we note here that if nib is neglected the resulting confidence intervals 
on Vr, gi are symmetric. The asymmetry between positive and negative couplings 
seen in Table Q] reflects the importance of the nib correction. It should also be noted 
that the nib dependence of the limits leads to a small systematic uncertainty, due to 
the uncertainty in nib. This is examined in detail in appendix lAl 



3 Angular asymmetries 

A simple and efficient method to extract information about the Wtb vertex is through 
angular asymmetries involving the angle 9* t between the charged lepton momentum (in 
the W boson rest frame) and the W + boson momentum (in the top quark rest frame). 
Alternatively, one may consider the angle 9a, between the charged lepton and b quark 
momenta in the W rest frame. Both approaches are equivalent since these two angles 
are related by 9\ + 9 a = 7r. (The determination of however, is simpler, because 
both momenta are measured in the same reference frame without any ambiguity in the 
boosts.) For any fixed z in the interval [—1, 1], one can define an asymmetry 

N(cos 9\ > z) - N (cos 9\ < z) 
z ~ N{cos9* e > z) + Njcos9j < z) ' ^' 

The most obvious choice is z = 0, giving the forward-backward (FB) asymmetry A FB 
[9,20]§| It is analogous to the FB asymmetries at LEP, which together with the ratios 
Rb, R c allow us to extract the couplings of the c and b quarks to the Z boson. The FB 
asymmetry is related to the W helicity fractions by 

A FB = ^[F R -F L }. (8) 

The measurement of this asymmetry alone is not enough to fully reconstruct the cos 9\ 
distribution. One can then think about other asymmetries for different values of z. 
The determination of Fi is easier if we construct asymmetries involving only Fr and 
F , or Fl and F . This is achieved choosing z = =f(2 2 / 3 — 1). Defining for convenience 
P = 2 1 / 3 - 1, we have 

2 = _( 2 2/3_i) ^ A t = A + = 3(3[F + (l + (3)F R ], 
z= (2 2 /3_ 1} ^ At = A- = -3/3[F + (l + l3)F L ]. (9) 

3 Notice the difference in sign with respect to the definitions in Refs. [9,20], where 9a is used. 
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Figure 4: Dependence of the asymmetries A + , and A FB on the couplings gL, gi 
and Vr, for the CP-conserving case. 



From both asymmetries and using F R + F L + F — 1, we obtain 

1 A^-PA + 
£r — - - + 



Fj 



L 



1-/3 3/3(1 — /3 2 ) 

1 A + - (5A_ 

1-/3 ~ 3/3(1 - /3 2 ) 



1 + /3 A + -A_ 

The three asymmetries A F b, A + , A_ are quite sensitive to anomalous Wtb inter- 
actions. Their SM values are A FB = -0.2225, A + = 0.5482, A_ = -0.8397, and 
their dependence on the non-standard couplings is shown in Fig. 01 Taking the ex- 
pected precision in their measurement from Refs. [15, 16] and assuming as central 
values the SM predictions, we obtain A FB ~ -0.223 ± 0.013, A + ~ 0.548 ± 0.010, 
A_ ~ —0.8397 ± 0.0033. Using e.g. the latter two, the helicity fractions can be deter- 
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mined as 



Fr = 0.0017 ±0.0071, 
F L = 0.2981 ±0.0167, 

F = 0.7002 ±0.0184. (11) 

The errors quoted take into account the correlation between the two measurements, 
which is determined writing the asymmetries in terms of the numbers of events in 
three bins: [-1, -(2 2 / 3 - 1)], [-(2 2 / 3 - 1), (2 2 / 3 - 1)] and [(2 2 / 3 - 1), 1]. We omit these 
details for brevity. The values extracted in this way are less precise than if obtained 
from a direct fit, but the method employed is much simpler too. The eventual limits 
which would be extracted from asymmetry measurements are collected in Table El A + 
exhibits the strongest dependence on gn and, if measured as precisely as it is expected, 
it would set the best limits on this coupling. On the other hand, A_ is the most 
sensitive to Vr and gz and sets the strongest bounds on them. The limits obtained 
from asymmetry measurements are competitive with those obtained from a direct fit 
to the cos 6} distribution. 

A+ A_ A FB 

V R [-0.15,0.15] [-0.056,0.11] [-0.12,0.15] 
g L [-0.12,0.082] [-0.057,0.026] [-0.092,0.062] 
g R [-0.019,0.018] [-0.024,0.022] [-0.027,0.025] 

Table 2: la bounds on anomalous couplings obtained from the measurement of angular 
asymmetries. 



4 Energy distributions 

The charged lepton energy in the W rest frame is fixed by the kinematics of the two- 
body decay W — > tv. Its energy in the top quark rest frame, denoted from now on by 
Ee, is related to the former by a Lorentz boost, and it is given by 

Et = ^(E w +\q\ cos 9* e ), (12) 

with \q\, given in Eq. (jHJ), the W boson momentum in the top rest frame and Ew 
its energy. Therefore, the angular distribution of the charged lepton in W rest frame 
determines its energy in the top rest frame. The maximum and minimum energies are 



10 



E 

Eqs. 



(E w + \q\)/2, E min = (E w — \q\)/2. The energy distribution is obtained from 



and (EG 
1 dT 

fdEp 



(E — E ■ ) 3 



[3(Ee — -Emin) 2 Fr + 3(-E max — Ft) 2 F L 



+Q(E max — E e ){Ei — E m - m ) F ] 



(13) 



The description of the top decay in terms of cos 9% or E e seems then equivalent, up to 
a change of variables. Any asymmetry built using cos 0|, defined around a fixed value 
z, can be translated into an equivalent asymmetry involving Ee, defined around a fixed 
energy E z = (E w + \q\z)/2, namely 



N{E £ > E z ) - N(Ee < E z 



(14) 



N(E e > E z ) + N(E t < E z ) 

However, in contrast to what was demonstrated in section [2] for the angular distribu- 
tions, finite width corrections have a non-negligible influence on Eg. This can be seen 
in Fig. El where we plot the energy distribution calculated analytically for t and W on 
shell and from a Monte Carlo calculation including finite width effects. The values of 
the asymmetries A + , A_ and Afb, calculated analytically and numerically (the latter 
for the cos 9\ and Et distributions) are shown in Table El One can notice the larger 
influence of finite width corrections for energy asymmetries. 




Figure 5: Differential distribution in Eq. ([TBI within the SM, calculated analytically 
and with a Monte Carlo simulation. 



The expected precision in energy asymmetries is worse than for the angular ones 
(given in the previous section) as it might be expected: Afb — —0.223 ± 0.024, A + ~ 
0.548 ± 0.013, A_ ~ —0.840 ± 0.016. Therefore, in principle their study does not seem 
to bring any improvement from the experimental side. 
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Analytical Angular Energy 



A + 0.5482 0.5492 0.5529 
A_ -0.8397 -0.8393 -0.8339 
^fb -0.2225 -0.2212 -0.2166 

Table 3: Values of the asymmetries A + , A_, A FB obtained from the analytical expres- 
sion (first column) and from the Monte Carlo simulation, the latter from the measure- 
ment of the distributions of cos^l (second column) and Ei (third column). 



5 Spin asymmetries 

Additional angular asymmetries can be built involving the top spin. Top quarks are 
produced unpolarised at the tree level in QCD interactions, and with a very small 
O(10~ 2 ) transverse polarisation at one loop. However, the t and i spins are strongly 
correlated, what allows the construction of angular asymmetries at the percent level. 
The top (anti)quark spins are not directly observable, but influence the angular dis- 
tribution of their decay products. For the decay t — > W + b — > £ + z/6, qq'b, the angular 
distributions of X = £ + ,u,q,q',W + ,b (which are called "spin analysers") in the top 
quark rest frame are given by 

i dr i, 

[1 + axcosOx) (15) 



T d cos &x 2 

with 6x the angle between the three-momentum of X (in the t rest frame) and the top 
spin direction. The constants «j are called "spin analysing power" of X and can range 
between —1 and 1. In the SM, o>i+ = — 1, a v — a q — —0.32, ayy+ = —a^ = 0.41 at 
the tree level [21] (q and q' are the up- and down-type quarks, respectively, resulting 
from the W decay). For the decay of a top antiquark the distributions are the same, 
with a x = —ax as long as CP is conserved in the decay. One-loop corrections modify 
these values to = 0.998, atqi = 0.93, a u = —0.33, a q = —0.31, acw+ = — a& = 0.39 
[22,23]. We point out that in the presence of non-vanishing Vr, Ql or gn couplings the 
numerical values of the constants a x are modified, but the functional form of Eq. (fT5l) 
is maintained. We have explicitly calculated them for a general CP-conserving Wtb 
vertex as written in Eq. (pQ) within the narrow width approximation. They can be 
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written as a x = cl X /o,q, with 

a = [Vl + V|] (1 + x 2 w - 2x 2 b ~ + AvA + 4) ~ V2x 2 w x b V L V R 
1-^-2^-^-^ + 4) 
- \2x 2 w x h - 6x w [V L g R + V R g L ] (l - 4k ~ 4) 



+ 6x w x b [V L g L + V R g R ] (l + x%- x 2 



= [Vl - V 2 ] (1 + x 2 w - 2x\ - 2x 4 w + x 2 w x\ + xt) - 12x 2 w x b V L V R 

1-^-2^-^-^ + 4) 

- \2x 2 w x b g L g R - 6x w [V L g R + V R g L ] (l - x 2 w - 4) 
+ 6 Wfc [Vitf £ - Vr^] (l + 4^ - 4) + 12x^(Fj - <?£) 
+ ^ log [_ 6 ^^ + 6 ^^2 (1 _ xl) + x^v^] , 

a,,, = [Vl - V|] (1 + x 2 w - 2x 2 - 2a& + x 2 w x\ + 4) + 124^ 6 V L V R 

1-^-2^-^-^ + 4) 

+ 12:r^:r 6 g L g R + Qx w [V L g R + V R g L \ (l - x 2 w - x\) 
+ 6x w x b [V L g L - V R g R ] (l + x 2 w - x 2 b ) - 12x 2 w (Vl - g\) 

+ ^ log [ 6x * w vl _ (1 _ xl) _ 12x 3 wXb V L g L ] , 

\q\ hj W - \q\ 

a b = -2^ { [V L 2 - ^] (1 - 2^ - 4) +2[ji- <&] (l-^ + 4 

+2x w [Vl^h - Vr# l ] + 6x w x b [V L g L - V R g R )} , (16) 

and aw = —a b . We have checked that our expressions are compatible with the first- 
order expansions in Refs. [24,25]. Working in the helicity basis and neglecting small 
spin interference effects, so that the cross section factorises into production times decay 
factors, the double angular distribution of the decay products X (from t) and X' (from 
t) can be written as [26] 

~1 Z~7 T~ = I^ 1 + Ca x a x ,cos9 x cos6 x ,) . (17) 

a a cos u x a cos u x , 4 

The angles 9 x , 9 X i are measured using as spin axis the parent top (anti) quark momen- 
tum in the ti CM system. The factor 

c = v(t R i R ) + a(t L i L ) - a(t R i L ) - a(t L i R ) 
~ cr(t R t R ) + a(t L i L ) + cr(t R t L ) + cr(t L t R ) 
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is the relative number of like helicity minus opposite helicitv tt pairs, and measures 
the spin correlation between the top quark and antiquarko We note that due to 
P invariance of the QCD interactions, cr{t^L) = cr{tijtR)t and by CP conservation 
c(£r?r) = ct(^l^l)- This is the reason why terms linear in cos6>x, cos9 x , are absent in 
Eq. (fl~7l) . In other words, terms linear in the cosines are present only if the top quarks 
are produced with a net polarisation in the helicity basis, what does not happen in 
pure QCD production. The actual value of C depends to some extent on the parton 
distribution functions (PDFs) used and the Q 2 scale at which they are evaluated. Using 
the CTEQ5L PDFs [28] and Q 2 = s the partonic CM energy, we find C = 0.310. At 
the one loop level, C = 0.326 ± 0.012 [23]. 

Using the spin analysers X, X' for the respective decays of t, i, one can define the 
asymmetries 

_ iV(cos 9 X cos 9 X , > 0) - iV(cos 9 X cos 9 X , < 0) 
xX> = iV(cos 6 X cos 9 X , > 0) + N(cos 9 X cos 0*, < 0) ' ^ ' 

whose theoretical value derived from Eq. (fT71) is 

A xx , = ~Ca x a x , . (20) 

If CP is conserved in the decay, for charge conjugate decay channels we have a X ia x = 
Oi X Oi X ^, so the asymmetries A x , x = A xx , are equivalent. Therefore, we can sum both 
channels and drop the superscripts indicating the charge, denoting the asymmetries by 
A u >, A U £i, etc. (CP- violating effects will be discussed in the next section). In semilep- 
tonic top decays we can select as spin analyser the charged lepton, which has the largest 
spin analysing power, or the neutrino, as proposed in Ref. [29]. In hadronic decays the 
jets corresponding to up- and down-type quarks are very difficult to distinguish, and 
one possibility is to use the least energetic jet in the top rest frame, which corresponds 
to the down-type quark 61% of the time, and has a spin analysing power acj = 0.49 at 
the tree level. An equivalent possibility is to choose the d jet by its angular distribution 
in the W~ rest frame [27]. In both hadronic and leptonic decays the b (b) quarks can 
be used as well. 

In the lepton + jets decay mode of the tt pair, tt — > tvbjjb we choose the two 
asymmetries Ay, A U j, for which we obtain the SM tree-level values Ay = —0.0376, 



4 Other conventions in the literature (e.g. Refs. [23,27]) denote by — C what in our case is the 
product C axctx'- We prefer to keep the notation in Refs. [1,26] and separate the contributions from 
the production (C) and the decay (ax, ot x >) since the former is sensitive to new physics in the tt 
production process while the latter are sensitive to non-standard Wtb interactions. This decomposi- 
tion is not possible if non-factorisable radiative corrections to the production and decay process are 
included. Anyway, these corrections are expected to be small. 
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A U j = 0.0120. With the precision expected for their measurement at LHC [16], the 
measurements A ej ~ -0.0376 ± 0.0058, A uj ~ 0.0120 ± 0.0056 are feasible. The 
dependence of these asymmetries on anomalous Wtb couplings is depicted in Fig. [6] (we 
remind the reader that the y axis scales are chosen so that the range approximately 
corresponds to two standard deviations around the theoretical SM value). These plots 
are obtained using Eqs. ([T6]) . (f20]) . We have checked, using high-statistics Monte Carlo 
simulations, that finite width effects are rather small, so that Eqs. (TT6l) . (l20l can be 
used to make accurate predictions for spin correlation asymmetries. In the dilepton 
channel tt — > Ivbl'vb we select the asymmetries Aw, A y £/, whose SM values are Au> = 
—0.0775, A V £i = 0.0247. The uncertainty in their measurement can be estimated from 
Refs. [16,19], yielding A u , ~ -0.0775 ±0.0060, A vtf ~ 0.0247 ±0.0087. Their variation 
when anomalous couplings are present is shown in Fig.[6l We also plot (in this case with 
arbitrary y axis scales) the asymmetries An,, Abb, which can be measured either in the 
semileptonic or dilepton channel. Their SM values are A tb = 0.0314, A b b = —0.0128, 
but the experimental sensitivity has been not estimated as yet. We expect that it may 
be of the order of 10% for A ib , and worse for A b b- 

The comparison of these plots with the ones in previous sections makes apparent 
that, given the experimental accuracies achievable in each case, spin correlation asym- 
metries are much less sensitive to non-standard Wtb couplings. This implies that, if 
no deviations are found in the measurement of the helicity ratios pr,l and angular 
asymmetries A±, spin-dependent asymmetries can be used to test tt spin correlations 
in the production, without contamination from possible new interactions in the decay. 
In particular, this is the case of Aw and Agj, whose relative accuracy is better, 7.7% 
and 15%, respectively. The determination of the correlation factor C in Eq. ( fl~8l) from 
these asymmetries would eventually give 

Au, - C = 0.310 ± 0.024 (exp) +°; 0043 (5V R ) (5g L ) ™ 4 6 (5g R ) , 

A ej - C = 0.310 ± 0.045 (exp) t°; 0068 (5V R ) ±™Z (Sg L ) t° Zt (Sg R ) . (21) 

The first error quoted corresponds to the experimental (systematic and statistic) un- 
certainty. The other ones are theoretical uncertainties obtained varying the anomalous 
couplings (one at a time). The confidence level (CL) corresponding to the intervals 
quoted is 68.3%. The numerical comparison of the different terms in Eqs. (|2T]) also 
shows that A^ and Aw are much less sensitive to non-standard top couplings than 
A +1 A- and p R ^- It must also be noted that, since all asymmetries depend on the 
production mechanism through the common factor C, their ratios do not (to leading 
order), and hence they are clean probes for anomalous couplings. The precision in 
the measurement of asymmetry ratios is still to be determined, but at any rate it is 
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Figure 6: Dependence of several spin correlation asymmetries on the couplings g^, ql 
and Vr, for the CP-conserving case. 



expected to be worse than for spin-independent observables discussed in the previous 
sections. 

It is also interesting to study the relative distribution of one spin analyser from the 
t quark and other from the i. Let fxx' be the angle between the three-momentum of 
X (in the t rest frame) and of X' (in the i rest frame). The angular distribution can 
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be written as [23] 



1 da 1 

= -{1 + Da x ax> cos ip X x>) , (22) 

a a cos (fxx 1 2 



with D a constant defined by this equality. In our simulations we obtain the tree-level 
value D = —0.217, while at one loop D = —0.238 [23], with a theoretical uncertainty 
of ~ 4%. Corresponding to these distributions, we can build the asymmetries 

~ N(cosipxx' > 0) - N(cosip XX ' < 0) 1 , . 

Ax *' = TTl ^ n x , AT/ — 7TT = -Da x a x > ■ (23) 

N(cosipxx' > 0) + N(cos(pxx' < 0) 2 

For charge conjugate decay channels the distributions can be summed, since ax>oi x = 
otxotx' provided CP is conserved in the decay. The dependence of these asymmetries 
A xx > on anomalous couplings is (within the production x decay factorisation approxi- 
mation) exactly the same as for the asymmetries A X x' defined above, and plots are not 
presented for brevity. Simulations are available for Ay and A U t , whose theoretical SM 
values are Ay = 0.0527, Aw = 0.1085. The experimental precision expected [16,19] 
is A ej ~ 0.0554 ± 0.0061, A u > ~ 0.1088 ± 0.0056. This is a better precision than for 
Ay and Ay, respectively, but still not competitive in the determination of the Wtb 
vertex structure]^ Instead, we can use them to test top spin correlations. From these 
asymmetries one can extract the value of D, obtaining 

A U , -> D = -0.217 ± 0.011 (exp)ir 031 (<%) {6g L ) (6g R ) , 

S3?, fe) ^axTs • 

(24) 



4 - D= -0.217 ±0.024 (exp) ^; 0047 Wr) ^xT« ( 5 9l) ™™ t {Sg 



The errors quoted correspond to the experimental (systematic + statistical) uncertainty 
and the variation when one of the anomalous couplings is allowed to be nonzero. As 
in the previous case, the measurement of a ratio of two asymmetries A XX ' provides a 
clean probe for anomalous couplings, but with a precision expected to be worse than 
for spin-independent observables. 



5 A special situation occurs if there is a fine-tuned cancellation between two nonzero Vr and gL 
couplings leading to small effects in W helicity fractions and related quantities. These cancellations 
are possible, and in such particular case the measurement of spin asymmetries like Aw and Au> (which 
are insensitive to gL but sensitive to Vr) or single top production may be used to obtain additional 
information about anomalous Wtb couplings. 
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6 Effect of complex phases in helicity fractions and 
spin asymmetries 



In the previous sections we have assumed that any non-standard Wtb couplings are 
real, either positive or negative. We have also pointed out that, if a non-zero coupling 
gn exists, its phase has an important influence on W helicity fractions and angular 
distributions determined by them. Complex phases in Vr, gi and gR influence the 
helicity fractions Fi through interference terms, which involve the real parts of these 
couplings (assuming Vl real). (Interference terms are the most important ones for 
small values of Vr, gi and gR, and for the latter coupling they are unsuppressed.) 
The maximum and minimum effects of anomalous couplings on Fi are obtained when 
they are real, negative or positive (not necessarily in this order). We show in Fig. [7] 
the values of the helicity fractions for fixed moduli and arbitrary phases of the new 
couplings, Vr = 0.1e l( ^ v R, g^ = 0.1 e 1 ^, gR = 0.1e"^« (one different from zero at a 
time), to illustrate the effect of the phases. The plot scales have been enlarged to cover 
all the range of variation of Fi, and the 2cr expected limits have been marked with a 
gray dashed line. 

For gL and Vr the deviations from the SM value are relatively stable under varia- 
tions of the phases, because the linear terms (which depend on the phase) and quadratic 
ones (which do not) are comparable in magnitude. Thus, the presence of a complex 
phase does not significantly affect the observability of the coupling. On the other hand, 
for gR the effect of the phase is dominant, and we notice that for phases 4> gR = ±n/2 
the helicity fractions are very close to their SM values, so that a purely imaginary cou- 
pling gR ~ O(0. 1) could remain unnoticed in an analysis of angular distributions. We 
also note that the plots are symmetric with respect to the y axis, because Fi depend 
on Re Vr, Reg^, RegR and the moduli. This also implies that complex couplings have 
the same effect on the helicity fractions in t and t decays. Therefore, the comparison 
in t and i decays of the angular distributions studied does not give any extra infor- 
mation regarding the complex phases, and further observables are needed in order to 
investigate this possibility. We have also analysed the phase dependence of the most 
interesting spin asymmetries, Ay, Ay, A(£> and Am (in semileptonic decays we consider 
W + — y W~ — > qq'). The effect of the phases is barely detectable, even with a 
greater experimental precision, and in any case the phase and modulus of an eventual 
anomalous coupling measured could not be disentangled. 
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Figure 7: Dependence of the helicity fractions F< on the phases of the anomalous 
couplings (see the text for details). 



7 CP-violating asymmetries 

In this section we examine whether the existence of complex phases in the anomalous 
couplings can be detected using CP-violating asymmetries. The possibility of a rela- 
tively large imaginary coupling is particularly intriguing, since angular distributions 
are very sensitive to this coupling provided its phase is not close to ±tt/2. For the other 
couplings, gi and Vr, the situation is not so dramatic, because the observability mainly 
depends on their moduli. 

The spin asymmetry [30] 

I rl cr{t R i R ) - a{t L t L ) 



(25) 



CP a{t R i R ) + a(t L t L ) ' 
is CP-violating, and vanishes at the tree level in QCD interactions. The top and antitop 
spins can be inferred using their decay products as spin analysers, in the same way as 
in section El We thus write 

RL N(cos6 x > 0,cos6>^, > 0) -iV(cos#x < 0,cos6^, < 0) 



A 



CP 



iV(cos X > 0, cos e x , > 0) + iV(cos 6 X < 0, cos X , < 0) 



(26) 
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Even with A^ vanishing in the production process, complex phases in the decay could 
in principle lead to an observable asymmetry. We have considered the dilepton channel, 
in which larger asymmetries are expected because of the higher spin analysing power 
of the charged leptons. (Other possibility to measure this asymmetry in the dilepton 
channel would be to consider the charged lepton energies [31].) We have found that, 
for anomalous couplings of order O(0. 1) and arbitrary phases, this asymmetry remains 
below the permille level, and with values consistent with zero within Monte Carlo 
uncertainty. Observation of such asymmetry would then unambiguously indicate CP- 
violating effects in ti production, which are possible, for instance, in two Higgs doublet 
models [30,32,33]. 

We also investigate triple-product asymmetries defined in the dilepton channel, 

Ti N(Tj > 0) - N(T t < 0) 

CP N(Ti > 0) + N(Ti < 0) ' 1 ' 

where the triple products T, are [32,34] 

Ti = e • (pe+ - Pi-) (pe+ xp t -)-e, 
T 2 = (p b - Pi) ■ (Pe+ t-Ptr), 

T 3 = (p t - pi) ■ {Pi + x p e -) . (28) 

The unit vector e is taken in the beam direction and the particle momenta follow ob- 
vious notation. Final state particle momenta can be measured in the laboratory frame 
or, if the kinematics of the event is completely reconstructed, in other reference system. 
For asymmetries built using Ti and T 3 we have found values O(10 -4 ), and compatible 
with zero, taking anomalous couplings of order O(0. 1) with arbitrary phases. On the 
other hand, we have found that AJ P is sensitive to a gn coupling of this magnitude. 
The asymmetry is larger if the charged lepton and b quark momenta are measured 
in the respective rest frames of the decaying top quarks. However, its observability 
will depend on systematic errors associated to the reconstruction, which have not been 
estimated as yet, and may be better when defined in the laboratory frame. The asym- 
metry A T q P in both reference systems as plotted in Fig. El for couplings Qr = 0.05, 
gn = 0.1 with arbitrary phases. 

Other asymmetries discussed in the literature are based on the differences 

Ai = — E(- , 

A 2 = Pi- pt+ ~ Pt ■ Pi- , 

A 3 = cos 9(+ — cos 6i- . (29) 

These quantities do not involve the product (pe+ x p t - ) , and then they can be measured 
in the semileptonic channel too [34]. The asymmetries are built as in Eq. f T27l) . and 
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Figure 8: Dependence of the CP asymmetry A^ p , defined in top quark rest frame (left) 
and laboratory system (right), on the phase of g R . 



take values below the permille level for anomalous couplings of the size considered in 
this section. For completeness, we have also considered the P-violating asymmetry 



I rl _ v(tRtL) ~ cr(t L t R ) 



Ap <j(t R t L ) + a(t L t R ) ' ^ 
which is measured using the charged leptons as spin analysers, 

aRL = n(cos e x > o, cos e x , < o) - iv(cos e x < o, cos e x , > o) , 

p N(cos9 x > 0,cos9 x , <0) + N(cos9 x <0,cos9 x , > 0)' [ ' 

We have found Ap L ~ 1.5 x 10~ 3 , and insensitive to anomalous couplings O(0. 1). This 
asymmetry can be sizeable in SM extensions [35]. 

We finally emphasise that, even in the cases where they are insensitive to complex 
anomalous Wtb couplings, the CP-violating asymmetries studied here are still very 
useful to disentangle CP violation in the production and the decay. If these asymmetries 
are found to be non-vanishing, they clearly signal CP violation in tt production. On 
the other hand, the imaginary part of Qr can be probed using A^ p . 



8 Summary 

New physics, if it exists close to the electroweak scale, may manifest itself through non- 
standard top interactions. In this paper we have discussed top pair decays at LHC as a 
probe of the Wtb vertex. We have examined angular and energy distributions, as well 
as asymmetries, involving or not the top quark polarisation. Among the observables 
discussed, the best sensitivity to anomalous Wtb couplings is given by the W helicity 
fractions Fj = Ti/T, i = L,R,0, and related observables. We have obtained analytical 
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expressions for Fi, for a general CP-violating Wtb vertex with the top quark and W 
boson on their mass shell, and keeping a non-zero bottom quark mass. We have shown, 
comparing with exact numerical results, the high accuracy of this approximation when 
studying angular distributions. We have also pointed out the importance of keeping 
the bottom mass in the calculations, in contrast with previous studies in the literature. 

W helicity fractions can be extracted from a fit to the charged lepton angular dis- 
tribution in W rest frame. The same analysis can be used to determine the helicity 
ratios p^i = Tr^/Tq, which are actually more sensitive to Vr and gi,-type anomalous 
couplings, given the experimental uncertainties (dominated by systematics already for 
a luminosity of 10 fb" 1 ) associated to each observable. A simpler method to probe the 
Wtb vertex, without the need of a fit to the charged lepton distribution, is through 
angular asymmetries. We have introduced two new asymmetries A + and A_, in addi- 
tion to the £W (or lb) forward-backward asymmetry A FB previously studied [9]. These 
new asymmetries allow us to: (i) obtain more precise bounds on anomalous couplings 
than Afb, comparable with those obtained from p^i and Fi, and even better for a Qr 
coupling; (ii) determine the helicity fractions with a fair accuracy without fitting the 
charged lepton distribution. The W helicity fractions determine the charged lepton en- 
ergy distribution in top rest frame as well. Energy asymmetries can be built, but they 
are less suited for the study of anomalous couplings because the approximation of con- 
sidering the top quark and W boson on shell is worse, and experimental uncertainties 
on energy asymmetries are larger. The best limits found, using single measurements, 
are 



Limits can be improved by combining the measurements of A± and Pr,l- The theoret- 
ical predictions for these and other observables have been implemented in a computer 
program TopFit, which allows to extract combined limits on anomalous couplings from 
a given set of observables, following the statistical approach outlined in appendix [Bj 
Detailed results including the correlation of the various observables (computed from 
Monte Carlo simulations) are beyond the scope of this paper, and have been presented 
elsewhere [15]. 

Spin correlations and spin-dependent asymmetries probe not only the Wtb inter- 
actions but also the dynamics of tt production. Their study is very interesting from a 
theoretical point of view, because they are sensitive to e.g. the exchange of a scalar 
particle in s channel [36] or anomalous gtt couplings [37]. It is then crucial to disentan- 



- 0.046 < g L < 0.013 

- 0.019 <g R < 0.018 



- 0.029 <V R < 0.099 



(pr) > 
(pr) , 



(32) 
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gle new physics in the production from possible anomalous Wtb couplings. This could 
be done, for instance, considering ratios of spin correlation asymmetries or, even bet- 
ter, using the strict bounds on anomalous couplings obtained from top decay angular 
distributions. 

The dependence of spin correlation asymmetries on Wtb anomalous couplings oc- 
curs through the "spin analysing power" constants of top quark decay products. We 
have calculated these constants for a general CP-conserving Wtb vertex, in the nar- 
row width approximation. It has been shown that the sensitivity of spin correlation 
asymmetries to top anomalous couplings is much weaker than for helicity fractions and 
related observables. Then, we have set explicit limits on the variation of two factors C, 
D (which measure the tt spin correlation) due to possible anomalous couplings not de- 
tected in other processes, i.e. within the ranges in Eqs. (l32l) . The possible variation in 
C, D is much smaller than the experimental precision expected for their measurement, 

C = 0.310 ± 0.024 (exp) 1°; 0043 (8V R ) (5g L ) ™ 4 6 (8g R ) , 

D = -0.217 ± 0.011 (exp) +_° 00i7 WO tSSS 6 - fe) t° 6 °7-s OM . (33) 

Hence, any deviation observed experimentally should correspond to new physics in the 
production. On the other hand, in ratios of two spin asymmetries A X x< [Axx>) the 
common factors C (D) cancel, and thus the ratios can cleanly probe non-standard top 
couplings. These observables have also been implemented in the computer program 
TopFit, and estimates for their expected precision will be presented elsewhere. 

Finally, we have addressed the possibility of complex anomalous Wtb couplings g^, 
g R ,V R . Complex phases in these terms influence helicity fractions and related quantities 
via the interference with the dominant SM coupling Vl (which we have normalised to 
unity). For V R and g^, quadratic and interference terms have the same magnitude, 
and the effect of phases is not very relevant. For g R , however, the interference term 
dominates, and the dependence on the phase is very strong. One finds that a g R 
coupling with a phase close to ±7r/2 has little effect on angular distributions, and even 
with a relatively large modulus it could remain unnoticed in such analyses. The same 
has been found for spin correlation asymmetries. However, we have shown that a CP 
asymmetry based on the triple product (pi, — pi) ■ (p£+ x p e -) is sensitive to a complex 
g R , taking values up to ±2% for g R = q=0.1i. If this asymmetry can be measured at 
LHC with a precision below the percent level, it could help to measure or bound g R . 
The remaining CP asymmetries analysed are very small, and insensitive to anomalous 
couplings of this size. Therefore, they can be used to isolate CP violating effects in 
tt production [32,33]. On the other hand, single top production at LHC can probe 
the Wtb interaction, and B or super-I? factories, with precise measurements of CP 
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asymmetries e.g. in b — > 57, might also give indirect evidence for (real or complex) 
anomalous Wtb couplings, helping to determine the structure of this vertex. 
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A Effect of rrib in the helicity fractions 

As it can be observed in Eqs. (j2J), interference terms involving Vr (or ql) and the 
dominant SM coupling Vl are proportional to Xb = mb/m t . These terms are of equal 
size as the quadratic terms for small Vr, g^, and cannot be neglected in the analysis. 
To illustrate their importance, we plot in Fig. [9] the dependence of the three helicity 
fractions on the anomalous couplings, for mj = 4.8 GeV and neglecting m&. The 
differences are apparent for Fr, and for F we have the extreme situation that the only 
dependence of this quantity on Vr is through the Xb term. 

The m,b dependence of the limits generates a small uncertainty due to the uncer- 
tainty in nib, for which we use the b quark pole mass. This b mass definition has an 
ambiguity of the order of Aq CD ~ 220 MeV (for other definitions the uncertainty is 
smaller). The variation of the limits in Eq. (|32l) when m& is taken as 4.8 GeV ±Aq C d 
is presented in Table [4] (we display additional digits in order to better illustrate the 
variation). The effect is larger for Vr and gi, as it is expected from the discussion in 
section El Nevertheless, the uncertainty only amounts to a few percent. 
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Figure 9: Dependence of the helicity fractions Fi on Vr and gz, for m b = 4.8 GeV and 
neglecting the b quark mass. 



Lower limit Upper limit 

Coupling m b - A QCD m b m b + A QCD m b - A QCD m b m b + A QCD 

V R -0.0303 -0.0293 -0.0281 0.0977 0.0994 0.1077 

9L -0.0444 -0.0456 -0.0461 0.0140 0.0135 0.0127 

g R -0.0192 -0.0194 -0.0192 0.0181 0.0180 0.0178 



Table 4: Influence of m b in the limits in Eqs. (1321 ): variation when the b quark mass is 
taken at its central value m b = 4.8 GeV, or adding and subtracting a small uncertainty 
A QCD = 220 MeV. 

B Extraction of limits from observables 



The derivation of limits on the anomalous couplings from the measurement of the 
experimental observables discussed has to be done with special care, due to the non- 
linear dependence of the latter on the former. In this appendix we explain the method 
we have used to obtain our limits. 
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Let us denote by O a generic observable, e.g. an angular asymmetry, and x an 
unknown parameter (in our case an anomalous coupling) upon which this observable 
depends, and for which we want to obtain a confidence interval. O is experimentally 
measured and is assumed to obey a Gaussian distribution (with mean and standard 
deviation given by its measurement). However, if the dependence 0(x) is non-linear in 
the region of interest, the probability density function (p.d.f.) derived for the parameter 
x will no longer be a Gaussian and a Monte Carlo method must be used to determine 
a confidence interval on x. 

We determine the p.d.f. of x numerically, using the acceptance-rejection method: we 
iteratively (i) generate a random value (with uniform probability) X{ within a suitable 
interval; (ii) evaluate the probability of 0(xi), given by the p.d.f. of O; (iii) generate an 
independent random number r j (with uniform probability) ; and (iv) accept the value 
if the probability of 0(xi) is larger than r,. The resulting set of values {x,} is distributed 
according to the p.d.f. of x given by the measurement of O. The determination of a 
central interval with a given confidence level (CL) 7 is done numerically, requiring: (a) 
that it contains a fraction 7 of the total number of values {xj}; (b) that is central, i.e. 
fractions (1 — j)/2 of the values generated are on each side of the interval. 

We have applied this method to obtain the limits on Tables CD and keeping only 
one of the couplings non-vanishing at a time. We point out that: 

1. The dependence on gn of the observables F i7 Pr,l, A± and A FB is approximately 
linear, as it can be observed in Figs. EHU Therefore, the limits on this coupling 
can be approximately obtained directly from these plots using the method in 
Refs. [16,19]: for a given observable O, intersecting the plot of O(gn) with the 
two horizontal lines O = O exp ± AO, which correspond to the la variation of O, 
gives the la interval (with a 68.3% CL) on g^. 

2. The dependence on g^ and Vl is highly non-linear (the region of interest is at the 
extreme of a quadratic function) , and appreciable differences are found between 
the Monte Carlo and the intersection methods. For example, the "lex" limit on 
Vr obtained with the intersection method from the (hypothetical) measurement 
p R ~ 0.0005 ± 0.0026 is -0.051 < V R < 0.12. However, this interval has a 
confidence level of 85.6%, and the true 68.3% central interval obtained from the 
same measurement with the Monte Carlo method outlined above is —0.029 < 
Vr < 0.099. Although overcoverage is not as bad as undercoverage, it is quite 
desirable that confidence intervals have exactly the CL they are supposed to have. 

A similar procedure is applied to estimate the theoretical uncertainties in Eqs. (|2"TT) 
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and pit) due to possible anomalous couplings. 
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